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Abstract 

A large deviation technique is used to calculate the microcanonical entropy function 
s(y,m) of the mean-field 99 4 -model as a function of the potential energy v and the 
magnetization m. As in the canonical ensemble, a continuous phase transition is 
found. An analytical expression is obtained for the critical energy v c (J) as a function 
of the coupling parameter J. 
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1 Introduction 



Statistical physics is the branch of theoretical physics which provides a mi- 
croscopic basis for the macroscopic theory of thermodynamics. Equilibrium 
statistical physics makes use of statistical ensembles, like the microcanonical 
or the canonical one. The basic thermodynamic function in the microcanonical 
ensemble is the entropy s(v, m, . . . ) as a function of some energy v, magneti- 
zation m, and/or further arguments, whereas in the canonical ensemble this 
role is played by the canonical free energy f(T, h, . . . ) as a function of the 
temperature T, a magnetic field h, and/or further arguments. The computa- 
tion of thermodynamic functions from statistical mechanics is a formidable, 
in most cases impossible task in general, but, as a rule of thumb, it is more 
difficult in the microcanonical ensemble. This is the reason why, although the 
microcanonical ensemble is the most fundamental one of the statistical ensem- 
bles, one often resorts to the calculation of the canonical free energy / instead. 
This is a reasonable way to proceed, as long as both options yield equivalent 
results. If this is not the case, one speaks of ensemble nonequivalence, and 
the occurrence of non-concave microcanonical entropy functions in systems 
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with long-range interactions is a hallmark of this situation. Then, unlike in 
standard thermodynamics where the concavity of the entropy is assumed, the 
entropy s cannot be obtained from the free energy / by means of a Legendre 
transform. 

Only in recent years the community of statistical physicists became aware 
that large deviation theory, a branch of probability theory, is a powerful tool 
to compute, among others, microcanonical quantities, at least for some classes 
of systems (see [1] for a mathematical treatise on large deviation theory with 
applications to statistical mechanics, or [2] for a less formal presentation of 
the subject). Explicit calculations of microcanonical quantities by means of 
large deviation theory can be found in [3] for some models of physical interest 
(the infinite range Potts model, the Hamiltonian mean-field model, and the 
Colson-Bonifacio model of a free electron laser). 

We start by briefly reviewing some results of large deviation theory and their 
relation to statistical mechanics in Sec. 2. In Sec. 3, the microcanonical entropy 
of the y9 4 -model without interaction is calculated using large deviation theory. 
From this result, the entropy of the mean-field <£> 4 -model in the presence of 
interactions is deduced. The physical implications of the result are discussed, 
followed by some conclusions in Sec. 4. 



2 Large deviation theory, a tool in statistical mechanics 



This section contains a brief, informal review of some statements of large devi- 
ation theory and how they can be applied to statistical mechanics. Consider a 
sequence of independent, identically distributed (i.i.d.) random variables G 
H d with probability distribution p(X), the empirical mean & N = Xl+ '^ +XiV , 
and the mean (or expectation) value \x := E[X] := / dX X p(X). The law 
of large numbers tells us that, in the limit of large N, the probability P of 
an event with (5 at 7^ A* converges to zero. The form of this convergence is the 
central issue of large deviation theory, thus allowing predictions about very un- 
likely events. Under some prerequisites on the random variables X i: Cramer 
[4] found that, in the limit of large N, the convergence is of exponentially 
decaying form in N, 

-I(x)= lim ^-lnP(e N e \x,x + dx}). (1) 

The so-called rate function I(x) depends on the value x the empirical mean 
has to meet. Comparing (1) to the definition of the Boltzmann entropy, 

s(x) = lim In / dX^ 5 \x - 6 N (x N )] , (2) 
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we see that s(x) is identical to —I(x) up to a physically irrelevant additive 
constant. Under suitable conditions on X*, large deviation theory gives an 
instruction how to derive the rate function (and therefore the micro canonical 
entropy): 

s(x) = -(t x ,x)+\nE[e^>% (3) 

where (, ) denotes the usual scalar product and t x (x) is obtained by solving 
the conditional equation 



xE[e 



E 



(4) 



3 Application to the mean-field <£> 4 -model 



The mean-field y9 4 -model we want to study is characterized by a potential 
energy function of the form 

t / N \ 2 N 

^) = -^(X>) +Ea^"|^), (5) 

where N is the number of particles, J is a coupling constant, and <p = 
(</?!, • • • , ip N ) denotes the position in configuration space. Our aim is to com- 
pute the microcanonical configurational entropy of this model. 



3.1 The tp A -model without interaction 



First, we look at the y9 4 -model without interaction, where N particles move 



independently in an on-site potential of the form \ipj — \tp\. The potential 



energy is 

N 



i=i 



The first step in calculating the configurational entropy is to translate our 
problem into the language of the previous section. Considering the coordinates 
<fi G R as random variables with flat distribution on some finite interval, 

else, 



other random variables 
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with potential energy and displacement of the i-th particle as components can 
be constructed from the </?j. This form of (8) is motivated by the fact that 
it allows to express the potential energy function (5) by the empirical mean 
& N (X) = (z N ,m N ), where 

M<p) = jrE\<pi-ltf, (9) 

1=1 

1 N 

™n{<p) = -rj^tPi- (10) 

JV i=l 

The associated macroscopical variables x = (z, m) are the mean potential en- 
ergy per particle z and the mean magnetization (or displacement) per particle 

TO. 



Then, following Eqs. (3) and (4), the entropy of the non-interacting <y? 4 -model 
can be written as 



s(z,m) = -t z z -t m m + \n d<p e tmlp+tz ( l ^ 4 ~^ 2 \ (11) 

J-fc 



where t z = t z (z, to) and t m = t m (z, to) are obtained by solving the conditional 
equations 



Z — \ uJi{t Z) t m ) — jLJ2(t z ,t m ), 
TO = Ui(t z ,t m ), 



(12) 



containing the integrals 



J^ c dip ip k e^v+tzi^ 4 -^ 2 ) 

Uk(tz,t m ) — — t : — — /i 4 i 2 \ ■ (13) 

f+£ c dip e' m ^ + *H4^ 4 -2^) 

We let p c tend to infinity in the following. 

Solving (12) for t z and t m seems not feasible analytically. However, besides 
a numerical solution (see Fig. 1), some general properties of s(z,m) may 
be of interest. From large deviation theory we know that s(z,m) is a con- 
cave, infinitely many times differentiable function. Its domain can be deter- 
mined by the following argument: The energy V(<pi) = \ <pf — \ <pf of any 
single particle is not bounded above, but it is bounded below by — |. There- 
fore also the mean potential energy z is restricted to the interval [— |,oo[. 
For any value of z ^ |, there exists a maximal and a minimal value of the 
magnetization m max {z) = —m min (z) which is attained when ipi = <pj for all 
So m max/min (z) is determined by z = \ m A max/min - \m 2 max/min . Starting 
from a microstate with (z, m max ), we can always find another microstate with 
to G [— m max , m max ] by changing the sign of an arbitrary fraction of the val- 
ues of the ipi. For an illustration of the domain of s(z,m) see Fig. 1. For any 
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Fig. 1. Entropy s(z,m) (left plot) and its domain (right plot, gray hatched area) 
of the (£ 4 -model without interaction from a numerical evaluation of (11). 

fixed value of z ^ t, the maximum of s with respect to m is located at zero 
magnetization. 



3.2 The ip 4 -model with mean-field interaction 



Turning to the study of the interacting model, we can now profit from our 
particular choice of the random variables Xi in Eq. (8), resulting in the macro- 
scopic variables (z,m). From the form (5) of the potential energy function V, 
the mean potential energy per particle can be written as 

v N ((p) = z N (ip) - j m N (Lp) 2 , (14) 

and for the corresponding macroscopic variable v the equality 

v = z — ^ m 2 (15) 

holds. As a consequence, the entropy s(v,m) of the mean-field y9 4 -model is 
obtained from the entropy of the non-interacting model by a simple transfor- 
mation of variables, 

s(v , m) = s(v + 7j m 2 , m). (16) 

A plot of the resulting function is shown in Fig. 2. For large enough fixed 
values of the potential energy v, the maximum of s with respect to m is 
again located at zero magnetization. Below a critical value v c of the potential 
energy, however, this is not the case anymore, indicating the occurrence of a 
phase transition. The transition is found to be a continuous one, and v c is 
determined by a change of curvature of s with respect to m along the line of 
zero magnetization, 

d 2 s(v c , m) 
dm 2 



= 0. (17) 

m=0 
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Fig. 2. Entropy s(v , m) of the mean-field (/? 4 -model for coupling constant J = 1 
(left plot) and its domain (right plot, gray hatched area). s(v,m) is obtained by a 
deformation (variable transformation) of the entropy s(z, m) of the non-interacting 
model (see Fig. 1). Below a critical value v c of the potential energy v, the maximum 
of s with respect to the magnetization m is located at a non-zero value of m. 

From the results of Sec. 3.1, the critical potential energy v c in dependence of 
the coupling constant J can be written as 



v c {J) 



where t Zc is denned implicitly by 



J t Zc uj 2 (t Zc ,0) = -1. 



;is) 



(19) 



It follows immediately that v c (l) = and, since t z < 0, that v c (J) > for all 
J > 0. The integral cl>2(£ 2c ,0) in (19) can be rewritten in terms of modified 
Bessel functions of the first kind Ik, yielding 
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+ h 



(20) 



This allows to derive asymptotic expansions of v c (J) for small and large pos- 
itive J, respectively, from the known expansions of these Bessel functions. 



Small J > 0: Inserting the asymptotic expansion of (20) for large negative 
t Zc into (19), solving for t Zc (J) and substituting into (18) finally yields 

v c (J) = -\ + U + 0(J 2 ). (21) 



Large J > 0: Analogously, an expansion of (19) for small negative t Zc yields 
an expression for the critical energy in the limit of large J, 

v c (J) = a* J* - (2 a 2 - I) J + (5f - § + ^) + O (l) , (22) 
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with a = r(|)/r(|) « 0.338. 

Together with the numerical result from the evaluation of Eqs. (18) and (19), 
these two asymptotic expansions are plotted in Fig. 3. 
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Fig. 3. Critical energy v c as a function of the coupling constant J from numerical 
computation (dots), confronted with the two asymptotic expansions (21) and (22). 



4 Conclusions 

The microcanonical entropy s(v, to) as a function of the potential energy v and 
the magnetization to of the mean-field y9 4 -model has been calculated using a 
large deviation technique. For such long-range interactions, and in contrast 
to standard thermodynamics, the entropy need not be a concave function. In 
fact, for a ferromagnetic coupling J > 0, a non-concavity in s is found, and 
therefore, when comparing this result to the canonical free energy f(T, h) as 
a function of the temperature T and the magnetic field h, nonequivalence of 
ensembles is observed. For any fixed value of v above a critical energy v c , the 
maximum of s(v, to) with respect to to is attained at to — 0, whereas below v c 
it is attained at m ^ 0, corresponding to a continuous order phase transition 
at v c . An exact, implicit expression for the critical energy v c (J) as well as 
expansions for small and large positive J are derived. 

Maximizing s(v,m) over m gives a concave entropy function 

s(v) — max [s(v, to)] , (23) 

showing ensemble equivalence when compared to the canonical free energy 
f(T) = maXft[/(T, h)]. In fact, the critical energy v c (J) from our microcanon- 
ical calculation corresponds to the critical temperature T C (J) of the canonical 
result reported in [5]. 
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